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ABSTRACT 

The  disturbance  velocity  fields  due  to  translational  and  rotational 
motions  of  an  ellipsoid  in  a  uniform  stream,  constant  vorticity  and  constant 
rate-of-strain,  required  in  fundamental  studies  of  behavior  of  suspensions, 
have  been  obtained  by  the  singularity  method.  These  solutions  extend  earlier 
solutions  for  prolate  spheroids.  Although  equivalent  solutions  were  obtained 
by  Oberbeck  (1876),  Edwardes  (1892)  and  Jeffery  (1922)  by  separation  of 
variable  in  ellipsoidal  coordinates,  the  singularity  solutions  are  far  more 
simple  in  form.  Other  significant  results  obtained  by  the  singularity  method 
include  the  exposition  of  the  unified  structure  shared  by  the  three  boundary 
value  problems  and  the  construction  of  new  forms  of  the  Faxen  laws  for 
ellipsoids  through  application  of  the  reciprocal  theorem.  The  disturbance 
solutions  and  Faxen  laws,  the  basis  for  Smoluchowski ' s  (1911)  method-of- 
ref lections  technique,  are  used  to  calculate  hydrodynamic  interactions  between 
two  or  more  arbitrarily  oriented  ellipsoids.  In  particular,  mobility  problems 
are  solved  directly  to  order  R~5,  where  R  is  the  centroid-to-centroid 
separation  between  the  ellipsoids .  v 
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SIGNIFICANCE  AND  EXPLANATION 


The  calculation  of  hydrodynamic  interactions  between  nonspherical 
particles  is  needed  for  the  understanding  and  control  of  many  natural  and 
manufacturing  processes,  for  instance,  those  involving  sedimentation, 
colloidal  stability  or  suspension  rheology.  While  single-particle  solutions 
for  disks,  needles  (and  in  general,  ellipsoids)  are  available,  such  solutions 
are  not  easily  generalized  to  multi-particle  problems,  necessary  for  examining 
particle-particle  interactions.  New  forms  are  presented  here  which  greatly 
facilitate  this  task.  A  model  sedimentation  problem  involving  two  ellipsoids 
of  revolution  is  used  to  illustrate  the  general  technique.  Future  work  will 
cover  other  applications  of  this  technique. 
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NOTATION 


resistance  tensor  for  spheroids, 
largest  semi-axis  of  ellipsoid, 
intermediate  semi-axis  of  ellipsoid, 
resistance  tensor  for  spheroids, 
smallest  semi-axis  of  ellipsoid, 
spheroid  orientation  vector, 
rate-of-strain  tensor, 
eccentricity  of  the  generating  ellipse, 
force  exerted  by  the  fluid  on  the  particle, 
density  function  in  singularity  distributions, 
gravitational  vector. 

resistance  tensor  for  spheroids,  rank  »  3. 
Oseen-Burgers  tensor. 

vector  operator  in  singularity  solution, 
resistance  function  for  spheroids,  rank  -  K. 
unit  vector  normal  to  the  surface, 
pressure. 

density  function  in  singularity  distributions, 
center  to  center  separation  between  two  ellipsoids, 
radial  coordinate  from  particle  center, 
symmetric  part  of  the  stress-dipole  (stresslet). 
torque  exerted  by  the  particle  on  the  fluid, 
anti-symmetric  part  of  the  stress-dipole. 

U  particle  translational  velocity. 
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v  fluid  velocity. 

X  resistance  function  for  spheroids, 

x  Cartesian  coordinate, 

x  position  vector, 

x'  point  on  the  fundamental  ellipse. 

Y  resistance  function  for  spheroids, 

y  Cartesian  coordinate. 

Z  resistance  function  for  spheroids, 

z  Cartesian  coordinate. 

Greek  Letters 

aQ,  a'Q,  Constants  in  Jeffery's  (1922)  solution. 

80,  Sq,  b-  Constants  in  Jeffery's  (1922)  solution. 

Yq.  Y',  Y''  Constants  in  Jeffery's  (1922)  solution. 

5  identity  tensor. 

alternating  tensor, 
polar  angle  for  spheroids, 
fluid  viscosity. 

ellipsoidal  coordinate  (p  constant  gives  ellipsoidal  surface), 
stress  tensor. 

X  ellipsoidal  harmonic. 

Q  Dirchlet  potential 

Q  angular  velocity  of  fluid. 

Li 


particle  angular  velocity. 


labels  for  particles. 

refers  to  the  fundamental  ellipse 

indices  used  in  the  Einstein  summation  convention. 

label  for  multipoles  in  the  singularity  solution. 


label  for  the  n-th  reflection. 


ambient  field. 


SINGULARITY  SOLUTIONS  FOR  ELLIPSOIDS  IN  LOW-REYNOLDS-NUMBER  FLOWS: 

WITH  APPLICATIONS  TO  THE  CALCULATION  OF 
HYDRODYNAMIC  INTERACTIONS  IN  SUSPENSIONS  OF  ELLIPSOIDS 

Sangtae  Kim* 

1 .  INTRODUCTION 

Suspensions  of  nonspherlcal  particles  exhibit  non-Newtonian  behavior  through 
the  interaction  between  the  flow  field  and  Brownian  motion  (Giesekus  (1962), 
Brenner  (1972),  Hinch  &  Leal  (1972).  However,  rigorous  derivation  of  the 
material  functions  to  date  have  been  restricted  to  the  dilute  limit,  partly 
because  of  the  lack  of  information  on  multi-particle  hydrodynamic 
interactions.  Existing  information  on  particle-particle  interactions  is 
limited  to  interactions  between  prolate  spheroids  in  certain  geometries  such 
as  large  particle-particle  separations  (Wakiya  1965)  or  special  configurations 
(Gluekman  et.  al.  1971;  Liao  &  Krueger  1980).  Hydrodynamic  interactions 
between  oblate  spheroids,  despite  widespread  occurence,  e.g.  the  disk-shaped 
kaolinite  minerals  In  clay/water  suspensions,  have  received  even  less 
attention. 

The  first  steps  towards  a  method-of-ref lections  solution  of 
multi-ellipsoid,  hydrodynamic-interaction  problems  are  presented  here.  Our 
primary  goal  is  the  solution  of  problems  where  the  rigid-body  motion  of  the 
particles  are  to  be  determined,  given  the  external  forces,  torques  and  the 
ambient  velocity  field.  Defined  by  Batchelor  (1976)  as  mobility  problems, 
these  problems  appear  most  frequently  in  the  hydrodynamic  interaction  terms  of 
the  rheological  theories  mentioned  above. 
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Our  faith  in  the  method-of-ref lections  approach  is  based  on  the 
experience  with  spherical  particles  where  it  is  known  that  one  can  solve  the 
mobility  problems  accurately  with  a  surprisingly  small  number  of  reflections 
(see  for  example,  Felderhof  (1977)  and  Jeffrey  &  Onishi  (19814)).  This 
conclusion  appears  to  hold  as  well  for  prolate  spheroids.  Kim  ( 1 98M t> )  has 
determined  the  sedimentation  velocities  of  two  arbitrarily  oriented  spheroids 
accurate  to  order  R  ^  where  R  is  the  centroid-to-centroid  separation,  using 
only  two  reflections  beyond  the  isolated-particle  solution. 

The  method  of  reflections  used  here  follows  Smoluchowski  (1911).  Readers 
who  are  not  familiar  with  the  details  of  this  technique  are  referred  to  the 
discussion  in  Happel  &  Brenner  (1973).  The  method  is  summarized  as  follows. 
The  disturbance  velocity  field  generated  by  a  test  particle  (call  this 
particle-a)  will  modify  the  velocity  field  seen  by  other  particles  (for 
example,  at  particle-3).  We  call  the  disturbance  velocity  generated  by  3  in 
response  to  the  disturbance  from  a  as  the  "reflected  field  at  3".  This 
process  can  be  continued  indefinitely,  with  each  reflected  fields  as  the 
incident  field  at  higher  order  reflections. 

The  disturbance  velocity  field  of  an  isolated  particle  can  be  considered 
a  reflection  with  the  ambient  velocity  field  as  the  incident  field.  We  call 
this  the  zero-th  reflection.  The  first  reflection  generates  reflected  fields 
at  each  particle  with  the  zero-th  reflection  fields  from  all  other  particles 
■.s  the  incident  fields.  For  an  M-particle  suspension,  the  N-th  reflection 
generates  M*(M-1)f'  reflected  fields  from  Mx(M-1)N  1  incident  fields. 

We  shall  represent  the  reflected  field  by  a  multipole  expansion  with  the 
multipcle  moments  related  to  the  incident  field  by  Faxen  laws  (Rallison  1978, 
Kim  1983).  Thi3  approach  is  direct  and  simple  but  its  success  hinges  upon  the 
availability  of  the  Faxen  laws  for  the  lower  order  moments. 
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Brenner  (19611)  has  shown,  using  the  Lorentz  (1907)  reciprocal  theorem, 
that  Faxen  laws  can  be  constructed  if  one  knows  the  stress  distribution  in  the 
conjugate  velocity  problem;  i.e.  the  solution  for  translating  ellipsoids, 
rotating  ellipsoids  and  ellipsoids  in  a  rate-of-strain  field  are  required  for 
the  Faxen  laws  for  the  force  (Brenner  1964),  torque  (Brenner  1964)  and 
stresslet  (Rallison  1978).  Furthermore,  if  the  conjugate  solution  is 
expressed  in  terms  of  the  fundamental  solution  of  the  Stokes  equation,  (also 
known  as  the  Stokeslet)  then  as  noted  by  Hinch  (1977),  Brenner’s  (1964) 
procedure  reduces  to  the  simple  statement  that  the  Faxen  law  for  the  moment 
has  the  same  functional  form  as  the  conjugate  velocity  field.  An  explicit 
statement  and  proof  is  given  in  Kira  (1985a).  Therefore,  the  bulk  of  the 
present  work  is  directed  towards  finding  such  singularity  solutions  for 
ellipsoids. 

The  singularity  method  has  been  used  by  Chwang  &  Wu  (1974,  1975)  to  solve 
exactly  the  translational,  rotational  and  rate-of-strain  problems  for  prolate 
sphenoids.  These  solutions  are  the  conjugates  for  the  Faxen  force,  torque  and 
stresslet  law.  In  their  introduction,  they  review  the  history  of  the 
singularity  method,  including  the  pioneering  works  of  Lorentz  (1892),  Oseen 
(1927)  and  Burgers  (1938).  Other  early  applications  of  the  singularity  method 
are  the  works  on  slender-body  theory  by  Hancock  (1953)  and  Tuck  (1964).  Less 
ha3  been  done  on  "thln-body"  theory  but  recent  results  are  available  for  thin 
oblate  bodies  of  revolution,  e.g.  Barshinger  &  Geer  (1984).  As  stated  by 
Chwang  h  Wu  ( 1 97 5 ) ,  "through  these  investigations,  the  relative  simplicity  and 
effectiveness  of  the  (singularity)  method  have  gradually  become  more 
recognized".  However,  the  primary  difficulty  is  that  a  priori,  one  does  not 
know  che  tyoe  of  singularities  end  their  distributions.  In  fact,  one 


4 


objective  of  the  work  of  Chwang  &  Wu  (197*1,  1975),  Chwang  (1975)  was  the 
accumulation  of  a  class  of  exact  solutions  by  the  singularity  method.  Section 
2  of  the  present  work  adds  the  general  ellipsoidal  shape  to  this  collection. 

The  organization  of  this  paper  is  as  follows.  In  Section  2,  the 
ellipsoidal  solutions  of  Oberbeck  (1876)  and  Jeffery  (1922)  are  re-expressed 
as  singularity  solutions,  i.e.  In  terms  of  the  fundamental  solution  of  the 
Stokes  equation.  (Edwardes'  0892)  work  is  contained  within  Jeffery's 
solution).  The  form  of  the  singularity  solution  is  surprisingly  simple.  In 
addition,  the  singularity  method  reveals  a  unified  structure  which  is  not 
apparent  in  the  traditional  expressions  in  ellipsoidal  harmonics.  This 
structure  suggests  new  forms  for  the  velocity  representations  for  nonspherical 
particles  in  capillaries  and  other  bounded  domains.  The  special  case  of 
oblate  and  prolate  ellipsoids  of  revolution,  including  a  complete  discussion 
of  the  resistance  tensors,  is  provided  in  the  appendix.  Section  3  is  a 
discussion  on  new  forms  of  the  Faxen  laws  for  the  force,  torque  and  stresslet 
on  an  ellipsoid  with  applications  to  the  method  of  reflections  for  two 
ellipsoids.  The  hydrodynamic  Interactions  between  two  sedimenting  oblate 
spheroids  have  been  determined  and  are  compared  with  the  results  obtained  in 
Kira  (1985b)  for  prolate  spheroids. 
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2.  THE  SINGULARITY  SOLUTION  FOR  ELLIPSOIDS 

In  this  section,  we  will  derive  the  singularity  solution  for  an  ellipsoid  in 
Stokes  flow.  We  will  start  by  describing  the  classical  boundary  value  problem 
followed  by  a  description  of  the  singularity  solution.  The  derivation  is 
outlined  in  the  appendix. 

2. 1  Description 

Consider  an  ellipsoid  with  semiaxes  of  lengths  a,  b  and  c,  with  a£b2c.  The 
ellipsoid  surface  satisfies: 


The  governing  equations  for  the  velocity,  v,  and  pressure,  p,  are  the  Stokes 
equations  for  low-Reynolds-number  flow, 

-Vp  +  pV2v  *  0,  (2.2) 

where  u  is  the  viscosity;  and  the  equation  of  continuity  for  incompressible 
flow, 

V-v  -  0.  (2.3) 

The  boundary  conditions  are: 

1)  On  the  ellipsoid  surface,  v  is  equal  to  the  particle's  rigid-body  mot:'  . 

V  “  U  +  <0*  X  .  >  ;  : 

2)  As  |  x  |  -->  '-0,  v  approaches  the  ambient  velocity,  i.e., 

v  — >  v=U+Q*x  +  E*x, 

whe~e  U  and  u>  are  the  particle  translational  anc  -otational  velocities  a-.:  h 


m 

Q  and  E  are  the  uniform  stream,  ambient  rotation  and  rate-of-strain.  From 
(2. 4b),  it  follows  that  Q*  is  one-half  of  7*v",  the  ambient  vorticity. 

The  solution  to  this  problem  can  be  expressed  using  a  distribution  of 
i ( x-*’ )/(8iry) ,  the  fundamental  solution  of  the  Stokes  equation.  I,  the 
Oseen-Burgers  tensor  given  by 


I.,  =■  —6  .  ,  +  — 3x.x.,  with 

lj  r  lj  r3  i  j* 


and  its  pressure  field. 


2  p 

pj  “ 


(2.5a) 

(2.5b) 


satisfy  the  Stokes  equation  with  point  forcing, 

-  j  +  “  ~  8mp6^6(x),  (2.6a) 

3rd  the  continuity  equation,  -r— ij  -  0  (2.6b) 

dX  ^ 

(see  Happel  &  Brenner  (1973)  Chapter  2). 

It  is  now  claimed  that  the  disturbance  velocity  field,  v-v",  can  be 
tten  as: 


x)-v”(x)  -  I  L(  n)  •  j  j  f(n)(x’,y’)ll  +  )  ^  ( 8irp )  dx’dy' ,  (2.7) 

E 


n»1 


i  th 


,  /  s  (2n-1)  2n-3 

Cfn)U'y)  •  ' 


(2.8a-e) 


EE 


,  ,  r.  X*  y  i1/2 

q(x,y)  *  [1  "  tt  ~  trJ 


aE 


,  2  2,1/2 
ap  *  (a  -c  )  , 


.  20/2 

b£  *  (b  -c  )  , 


(n) 


(  -  F  if  n-1, 

\  5*V  +  X-7  if  n=2. 


following  is  a  description  of  the  terms  which  appears  in  the  solution. 
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:  the  orientation  trajectories  are  followed  from  9=0  (horizontally 

oriented  oblate  spheroids  and  vertically  oriented  prolate  spheroids)  the 
■orves  in  Figures  2  and  3  fall  into  two  groups,  depending  on  the  initial  value 
of  P .  For  both  oblate  and  prolate  spheroids,  if  R  (at  9=0)  exceeds  a  critical 
value,  then  the  particles  monotonically  drift  apart.  Their  orientations 
approach  asymptotically  a  limiting  value  of  9,  since  at  large  separations  u 
gees  to  zero.  However,  for  initial  values  of  R  less  than  the  critical  value, 
the  rotational  motion  is  sufficiently  large  to  cause  the  particles  to  rotate 
oevond  9=tt/2,  whereafter,  the  particles  drift  back  towards  each  other  along 
trajectories  which  are  mirror  images  of  the  outward  trajectories.  The 
sepa^atrix  which  starts  at  the  critical  value  of  R  has  the  asymptote  9"-it/2 
(horizontal  orientation). 

Ac  large  values  of  R,  the  trajectories  can  be  approximated  accuarately  by 
taking  just  the  leading  terms  on  the  right-hand-side  of  evolution  equations 
;.12)  and  (3‘i3).  These  approximate  equations  have  exact  solutions, 

oblate  spheroids :  1  _  F  “  ( 2/3)  ( 1  /YA  -  1/XA)(cos29  -  cos20j, 

R  Rq  0 

on 1  prolate  spheroids:  F  "  F  *  ( 2/ 3 ) ( 1 /XA  -  1/YA)(cos29  -  cos28  ). 

rt0  0 

A  A 

T r, resistance  functions  X  and  Y  are  defined  in  the  appendix. 

The  influence  of  the  aspect  ratio  is  seen  by  comparing  Figures  2  and  3 
(  o’’  aspect  ratios  of  10  and  2  respectively.  As  the  aspect  ratio  is  reduced, 

‘.he-  ~  eg  ion  occupied  by  periodic  trajectories  enlargens  and  the  trajectories 
■ t-a : ‘  “ti  into  the  vertical  lines  of  the  spherical  case.  Finally,  at  a  fixed 
;  'ti-~,  at  -  -.nger  hyd-c dynamic  interactions  between  oblate  spheroids 
'"‘•‘it  i -i  far  greater  reg'on  of  periodic  trajectories. 


Figure  3.  Evolution  of  orientation  and  separation  between  two  spheroids,  aspect 
ratio  =  2,  settling  as  in  Figure  1.  The  solid  and  dashed  trajectories 
are  for  oblate  and  prolate  spheroids  respectively.  The  separatrix  goes 
through  the  x  on  the  right  border. 
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.  now  consider  two  Inclined  oblate  spheroids  settling  with  their  axes 

lying  in  a  common  vertical  plane  (Figure  1).  At  all  times,  the  geometry  is 

sped* led  by  the  dimensionless  center-to-center  separation,  R/a  and  9,  the 
U 

polar  angle  between  d^  and  the  x-axis.  At  all  but  small  separations,  the 
two-  reflection  solution  provides  accurate  answers.  The  convergence  behavior 
is  similar  to  that  reported  in  Kim  (1985b)  for  prolate  spheroids. 

Tne  evolution  of  the  geometry  is  caused  by  the  anisotropy  in  the  mobility 
tensors  and  the  rotation  of  the  spheroids  about  the  y-axis  (which  for  oblate 
spheroid  is  co- incident  with  a  major  principal  axis).  Since  the  mobility  is 
less  in  the  axial  than  in  the  transverse  direction,  an  inclined  spheroid 
drifts  horizontally  as  it  settles.  At  the  same  time,  the  spheroid  rotation 
changes  the  orientation  of  the  axis.  These  two  effects,  under  the 
quasi-steady  assumption,  are  governed  by  the  dimensionless  equations  (with  R/a 
rewritten  now  as  R) 

• 

6  -  u  (R,e)  (3. '12) 

r  -  -2Ux(R,e)  (3.13) 

Figures  2  and  3  show  the  evolution  of  R  and  9  as  determined  by 
integrating  (3*12)  and  (3-13)  with  a  fourth  order  Runge-Kutta  routine.  The 
solid  lines  are  the  new  results  for  oblate  spheroids  and  the  dashed  lines  are 
the  earlier  results  for  prolate  spheroids.  The  plots  include  the  curve 

R  -  2(1  -  e*cos*9)1/2, 

for  contact  between  the  two  spheroids. 

*4.  This  ft  differs  from  the  one  used  in  Kim  (1985b)  for  prolate  spheroids  by 
-/2  in  order  that  the  same  value  of  9  in  the  two  problems  yields  identical 
cross-sections  in  the  x-z  plane. 
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The  contributions  to  the  sedimentation  and  angular  velocities  on 

ellipsoid  1  at  the  second  reflection  are  obtained  by  using  as  the  incident 

field  in  the  Faxen  laws.  The  method  of  reflections  result  for  is  now 
-5  -6 

accurate  to  0(R  ).  An  error  of  0(R  )  comes  from  the  neglected  quadrupole 

fields  in  v,j2  (for  which  the  Faxen  laws  are  as  yet  unavailable).  Thus  the 
solution  for  ellipsoids  has  been  developed  to  the  same  level  as  that  presented 
for  prolate  spheroids  in  Kim  (1985b). 

It  should  be  clear  from  the  steps  used  at  the  first  and  second 
reflections  that  in  general,  at  higher  order  reflections,  the  contribution  to 
the  sedimentation  velocity  of  ellipsoid  a  from  the  n-th  multipole  from 
ellipsoid  B  is  of  the  form: 


E  E 
a  B 


I'  *  (y&f  *  I<x-j!)/<8.„>  8 A  dA„, 


l2Va  $  8 


a  B 


where  L.  is  the  appropriate  raultipole  moment  on  ellipsoid  B  obtained  at  the 
previous  reflection.  An  analogous  procedure  can  be  followed  to  determine  the 
contribution  to  the  angular  velocity. 

In  actual  computatation,  the  integrals  over  the  fundamental  ellipses  were 
parametrized  with  the  elliptic  coordinates: 


x  -  aEpcos<j>,  y  »  b^psim}*. 

Three-point  Gaussian  quadratures  were  adequate  for  the  p-integration  (the 

~2  2 

Gaussian  quadratures  were  performed  after  the  change  of  variable,  p  »  1-p  ). 


Simpson's  rule  was  used  for  the  ^-integration. 


We  proceed  to  solve  the  sedimentation  problem  for  two  ellipsoids  by  the  method 
of  reflections  as  an  application  of  the  general  results  of  Section  2.  The 
zero-th  order  solution  at  particle-a  is  simply  the  velocity  field  generated  by 
an  isolated  ellipsoid  subject  to  an  external  force  F  .  For  example,  at 

-a 

ellipsoid  2, 


D^0)  -  (pA)_1 «F2  ,  (3.10a) 

-2  "  -a'JJ  f(1)<x,,y,)l1  ♦  ^c2q2?*}  I(x-x’)/(8iry)  dx'dy’,  (3.10b) 

E 

(Each  ellipsoidal  particle  in  the  suspension  has  its  own  fundamental  ellipse 
and  constants.  This  dependence  is  not  expressed  in  order  to  simplify  the 
notation).  The  contributions  to  the  sedimentation  and  angular  velocities  on 
ellipsoid  1  from  the  first  reflection,  and  are  obtained  by  using 

y2(x)  as  the  incident  field  in  the  Faxen  laws  for  ,he  translational  and 
rotational  velocities  on  a  force-free  and  torque-free  ellipsoid.  The  leading 
term  in  the  reflected  field,  «  is  3  Stokes-dipole  field, 

y21  (  x)  -  (S^.V).jj  f(2)(x’,y,){l  ♦  -gC2q2V2 }  I(x-x'  )/(8nu)  dx'dy', 

E 

with  the  stresslet  determined  from  the  appropriate  Faxen  law,  equation  (3  ','1 
The  first  reflection  at  ellipsoid  2  follows  in  a  similar  fashion,  and  the 
expression  for  the  analogous  dipole  field,  y12,  is  obtained  by  switching  ihe 


particle  labels. 


(3.5) 


T  -  pC  *  j  j  f(2)(x',y')  [  ^Vxv’Cx’ )  -  u  ]  dx'dy’ 

E 

+  yH: |j  f(2)(x’,y'){l  +-gc2q272}  e’(x')  dx'dy', 

E 

S1J  “  ^ijkjj  f(2)(x’*y,)^  ekl(5,)  dx’dy’  (3<6) 

E 

*  j II  f(2)(x,*y,)  [  jV-v"(x’)  -  0)  ]k  dx’dy’. 

E 

For  force-free  and  torque-free  particles,  these  results  can  be  rearranged 
into  the  following  expressions  for  the  translational  velocity,  rotational 
velocity  and  stresslet. 


|J  -  |j  f(i)(x'»y’){l  +  -^c*q2V2}  v“(x')  dx'dy', 

E 

oij  -  ||  f(2)(x',y')  CjVxy“(x,)]1  dx'dy' 

E 

*  <c")uHktj{j  *  Wv'l  v<*’> 

E 

S. <  -  m[m« ..  -  -  H  . ,(C_1 )  H  ,  „ ] 
ij  L  ijk2.  raij  ran  nklJ 


dx'dy' 


(3.7) 

(3.8) 


(3.9) 


*  ||  f(2)(x'»y,)l1  +  ^c*q2V2 }  ekjl(x')  dx'dy'. 

E 

Equations  (3.4)  to  (3.9)  reduce  to  the  appropriate  Faxen  laws  for  prolate 
spheroids  derived  by  Kim  (1985a)  in  the  limit  as  b  — >  c. 

Thus  the  force,  torque  and  stresslet  on  an  ellipsoidal  particle  in  and 
ambient  flow  field  v  are  obtained  by  integrating  the  ambient  velocity, 
vorticity  and  rate-of-strain  (respectively)  over  the  fundamental  ellipse, 
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j .  r  /tXEN  LAWS  FOR  ELLIPSOIDAL  PARTICLES 

A  correspondence  between  singularity  solutions  and  Faxen  laws  follows  as  a 
corollary  of  the  Lorentz  (1907)  reciprocal  theorem  (Brenner  1964,  Kim  1985a). 
The  new  forms  of  the  Faxen  laws  obtained  in  this  manner  are  more  useful  than 
earlier  infinite  series  expansions  derived  by  Brenner  (1964)  and  Rallison 
(1978)  when  the  higher  order  derivatives  of  the  velocity  field  are  not 
available. 

The  linear  relations  between  the  drag,  torque  and  stresslet  on  the 
ellipsoid  and  the  ambient  field  can  be  expressed  as: 


Ti  ■  “cu<n  '"’j  *  "HukEjk 


sij  '  “"ljki'kl  *  “Hkij<!1  ‘“>k 


(3.1) 

(3.2) 

(3.3) 


where  A,  B,  C,  H  and  M  are  material  tensors  whose  components  may  be  deduced 
from  equations  (2.9),  (2.10)  and  (2.11).  The  Faxen  relations  are 
generalizations  of  (3.1).  (3.2)  and  (3*3)  since  they  give  F,  T  and  S  in  any 
ambient  velocity  field  that  satisfies  the  Stokes  equations  over  the  unbounded 
domain.  We  apply  the  reciprocal  theorem  to  the  singularity  solutions  as  shown 
in  Kim  (1985a)  to  obtain  the  following  forms  of  the  Faxen  laws  for  the  force, 
torque  and  stresslet: 


F  «  yA-jj  f(1)(x’,y,){l  ♦  |c2q2V2}  v"(x')  dx'dy’  -  pA-U  , 


(3.4) 
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In  summary,  the  basic  results  are: 

1)  The  disturbance  velocity  field  for  a  translating  ellipsoid  (or  a  fixed 
ellipsoid  in  a  uniform  stream)  is  generated  by  a  distribution  of  stokes- 
lets  and  potential  doublets  over  the  fundamental  ellipse. 

2)  The  disturbance  fields  for  a  rotating  ellipsoid  (or  a  fixed  ellipsoid  in 
a  constant  vorticity  field)  and  for  a  stationary  ellipsoid  in  a  rate-of- 
strain  field  are  generated  by  a  distribution  of  rotlets,  stresslets  and 
Stokes-octupoles  over  the  fundamental  ellipse. 

3)  For  prolate  spheroids,  the  fundamental  ellipse  degenerates  into  a  line 
segment  from  one  focal  point  to  the  other  and  the  singularity  solutions  of 
Chwang  &  Wu  are  recovered.  For  oblate  spheroids,  the  fundamental  ellipse 
degenerates  into  a  circular  disk  with  a  diameter  equal  to  the  focal  length 
of  the  ellipse  of  rotation. 

In  all  cases,  the  density  functions  for  the  dominant  singularities  are  similar 

to  those  which  appear  in  analogous  problems  in  potential  theory. 


obtained  by  evaluating  the  following  harmonic  functions  at  X«0. 


,00 

x(X)  -  abc  [P(X)]“1dX,  (2.12) 

JX 

/CO 

a(X )  -  abc  [(a2+X)P(X)]_1dX,  (2.13) 

JX 

with  P(X)  -  [(a2+X)(b2+X)(c2+X)]1/2 

The  lower  limit  of  the  definite  integral,  X(x,y,z),  is  the  positive  root  of 

*1.  +  _£!_  +  _|L_  _ 

Pn  P+X  P^x 

The  functions  B(X)  and  Y(X)  are  obtained  by  successive  cycling  of  the 
dependence  on  a,  b  and  c.  The  '  functions  are  defined  by: 

o'(X)  -  (Y-B)/(b2-c2),  (2.1*0 

with  B'(X)  and  Y'(X)  defined  by  successive  cycling  of  the  dependence  on  a,  b 

and  c  (and  therefore  also,  a,  B  and  Y).  The  n  functions  are  defined  by: 

a"(X)  -  (b2B-c2Y)/(b2-c2),  (2. 15' 

with  B" ( X )  and  Yn(X)  defined  by  successive  cycling  of  the  dependence  on  a,  b 

and  c  (and  therefore  also,  a,  B  and  Y).  This  completes  the  description  of  the 

singularity  solution.  For  the  special  case  of  ellipsoids  of  revolution,  these 
constants  are  given  in  the  appendix  along  with  a  complete  description  of  the 
resistance  tensors.  Asymptotic  expressions  for  the  spheroidal  resistance 
functions  are  also  provided  for  slender,  flat  and  near-sphere  limits. 


3.  a,  B  and  Y  are  as  defined  in  Happel  &  Brenner  (1973)  and  differ  from 
Jeffery's  definition  by  a  factor  of  (abc).  This  also  holds  for  the  ’ 
"  functions. 
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an  elliptical  disk,  as  can  be  seen  by  looking  at  the  limit  c  -  0  in  equation 
(2.7). 

The  function  q(x,y)  which  appears  in  plays  a  prominent  role  in  the 

potential  theory  for  ellipsoidal  particles  (see  Miloh  (197*0).  In  fact,  in 

potential  theory,  q  1  is  the  requisite  charge  distribution  over  the 

fundamental  ellipse  which  generates  ellipsoidal  equipotential  surfaces. 

Chwang  &  Wu  (1975)  have  noted  that  the  distribution  of  Stokes  multipoles  in 

low-Reynolds-number  problems  is  similar  to  the  distribution  of  multipoles  in 

analogous  problems  in  potential  theory,  except  for  the  presence  of  additional 

2 

degenerate  raultipoles  (the  V  I  term)  in  equation  (2.7).  The  presence  of  such 
quadrupoles  (or  potential  doublet)  when  n-1  and  octupoles  when  n-2  in  (2.7) 
are  consistent  with  (and  in  fact  extend)  the  rules  stated  by  Chwang  &  Wu 
(1975)  for  prolate  spheroids. 

To  complete  the  solution,  we  must  relate  F,  T  and  S  in  terms  of  the 
knowns,  u“-U,  Q-u  and  E.  These  relations  are  found  in  Oberbeck  (1876)  and 
Jeffery  (1922),  some  of  which  are  shown  below.  Expressions  for  other 
components  can  be  obtained  by  the  well  known  mnemonic  of  cycling  the 
subscripts  x,  y  and  z,  and  the  dependence  on  a,  b  and  c. 

Fx  =  I6iruabc(x0  ♦  c^aV  (if-U^  (2.9) 

Tx  =  ^rMabc(b2e0  ♦  c2Y0)_1  [(b2  +  c2)(n”-u,x)  ♦  (b2-c2)l(Eyz  +  Ezy)  ]  (2.10) 

sxx '  vuab0  (2,Se„-«»Vo!k)"5,5  *  Yo°o  *  w" 

3  =3  =  ^truabc( a2an  ♦  b28_.)  ^  (2.11b) 

xy  yx  3  0  0 

”  [<a2-»2><VV  *  (vV[Tjr'  i(Exy.Eyx)]. 

Here,  Xq»  ^q.  Sq,  Yq,  a^,  8^,  Y^,  a^,  8^  and  Y”  are  constants  which  are 
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nrst  consider  l/n\  the  vector  operator.  For  n-1 ,  l/^*I  -  -F*I.  the 
Stokes  monopole  field,  where  F  is  the  force  exerted  on  the  ellipsoid  by  the 
fluid,  i.e. 

o *n  dA 

(2) 

For  n-2,  L  *1  -  [(S+T)*V]*I  is  the  Stokes  dipole  field,  where  S  and  T  are 
the  symmetric  and  anti-symmetric  parts  of  the  stress-dipole  on  the  ellipsoid, 

( o  *n)x  dA 

at  —  — 

S  is  called  the  stresslet  in  Batchelor  &  Green  ( 1 972 ) 1 .  T  is  related  to  the 
torque  exerted  on  the  ellipsoid  by  the  fluid,  by  the  usual  relation  between 
anti-symmetric  dyadics  and  pseudo-vectors, 

TiJ  “  "  2® iJkTk‘ 

the  interior  of  the  fundamental 


elliptical  disk  in  a  family  of 

semi-axes  of  the  fundamental 
2 

(2.8c)  .  The  density  function 
surface  singularity  distribution  for 


E(x’,y'),  the  integration  domain,  is 
ellipse, 


The  fundamental  ellipse  is  the  degenerate 
confocal  ellipsoids.  The  major  and  minor 
ellipse,  a£  and  b£,  are  given  in  equation 
f(n)Cx*,y*)  in  E(x',y’)  is  physically  the 


1 .  The  isotropic  part  of  the  symmetric  stress-dipole  usually  has  no  physical 
significance.  Batchelor  and  Green  (1972)  remove  this  degree  of  freedom  by 
setting  the  trace  of  the  stresslet  to  zero. 

2.  Hobson  (1955)  and  Miloh  (1974)  use  k  and  (k*-ha)^^  in  place  of  our  a_ 

and  b„. 

E 


-W_\' 
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Appendix  1 .  DERIVATION  OF  THE  SINGULARITY  SOLUTION 

We  show  here  that  the  singularity  solution,  equation  (2.7),  is  equivalent  to 
the  solutions  obtained  by  Oberbeck  (1876)  and  Jeffery  (1922),  i.e.,  the 
solutions  obtained  by  separation  of  variables  in  ellipsoidal  coordinates. 

The  proof  is  simple  once  certain  integral  representations  for  \{\)  and 
ft(A),  the  Dirchlet  gravitational  potential  for  a  solid  ellipsoid,  are 
established.  The  Dirchlet  potential  is  defined  bya: 

°<»>  ■  'ab«  ^  (i£r*  ’if-  <*•’> 

The  lower  limit  of  the  definite  integral,  A(x,y,z),  is  as  defined  earlier. 

The  required  Integral  representations  are: 

X  -  2abc  ||  f(1)(x',y')  jx-_xT|'  dx'dy’,  (A. 2) 

E 

ff  2  1 

n  -  -Ihrabc  jj  f(1)(x\y»)  q*(xV)  dx'dy’ ,  (A. 3) 

E 

with  f^j  given  by  equation  (2.8a). 

The  integral  representation  for  x  is  derived  in  Miloh  (197*0  in  the  more 
general  setting  of  representation  theorems  for  external  Lame  functions^,  F™. 
The  harmonic  x  is  related  to  the  Lame  functions  by  x(i)  "  2abcF°(p)  with 
p2~a2*X.  The  Integral  representation  for  ft  follows  from  the  representation 
for  x-  From  Kellogg  (1953). 

ft  -  x(i;u)/u  du.  (A. 4) 

The  parameter  u  is  introduced  by  replacing  a,  b  and  c  in  the  expression  for  x 


a.  This  definition  is  the  same  as  in  Happel  &  Brenner  (1973)  but  differs  from 
Jeffery's  by  a  factor  of  wabc. 

b.  These  functions  are  defined  in  the  extensive  treatise  by  Hobson  (1955). 


with  ua,  ub  and  uc.  According  to  (A. 4),  n  is  a  superposition  of  x  functions 
for  a  family  of  ellipsoids  imbedded  inside  the  original  ellipsoid.  This 
representation  is  one  way  of  demonstrating  that  a  is  a  harmonic.  The  desired 
result,  (A. 3),  can  be  obtained  by  inserting  (A. 2)  into  (A. 4)  and  performing 
the  u-integration  first. 

We  are  now  in  a  position  to  recover  Oberbeck’s  solution  as  given  by 
equation  (5*11.8)  in  Happel  and  Brenner  (1973)  for  an  ellipsoid  in  a  uniform 
stream  (streaming  in  the  x-direction  with  velocity  U*): 


~  2^  (x0  -0 


+cua2)  1  Irr  +  U*(  x,,+a,,a' 


V  fxll  -  v)  +  u' 
0  “0“  ;  ^X3x  XJ  U 


a1.  »,  2.-1 

vy  •  •  2^  <  VV  > 


3x 
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(A. 5) 


fV 

2ir 


,  2,-1 
(x0+V  5 


a*n 

3x3z 


+  U"(x0+a0a2)_1  xl£ 


The  distribution  of  Stokeslets,  I^tx-x’),  in  the  singularity  solution 

can  be  decomposed  as 


Iij(x-;*)  -  9Xi  -  -p^ij)  +  xj  3x ^  ' 


with  r 


|x-x' 


(A. 6) 


The  x_term  in  Oberbeck's  solution  is  obtained  by  taking  the  first  term  on  the 
right-hand-side  of  (A. 6)  and  recognizing  that  the  integral  over  the 
fundamental  ellipse  in  (2.7)  is  precisely  the  integral  representation  for  x* 
equation  (A. 2).  It  is  not  difficult  to  show  that  the  integral  of  the 
remaining  term  on  the  right-hand-side  of  (A. 6)  over  the  fundamental  ellipse  is 
related  to  the  second  derivatives  of  ft,  i.e., 


ft 

Ji 


f(1)(x\y’)  x J 
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3x, 


x-x’ 


32n 

3xi 3Xj  ' 


^rabc 


dx  ?dy’ , 


(A. 7) 
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with  the  temporary  notation,  a^-a,  a2*»b  and  a^»c.  This  completes  the 
transformation  of  the  Stokeslet  distribution. 

The  potential  doublet  satisfies 

V2I  -  -2VV(r-1) 

for  r*0.  Therefore,  the  integral  of  the  potential  doublet  over  the 

fundamental  ellipse  in  (2.7)  is  just  the  second  derivative  of  the  integral 

representation  for  n.  Thus  the  n-terms  in  Oberbeck's  solution  are  obtained  by 

combining  (A. 7)  and  the  distribution  of  potential  doublets  in  (2.7). 

By  a  similar  albeit  more  tedious  procedure,  one  can  relate  the  n*2  case 

in  the  singularity  solution  to  Jeffery's  (1922)  solution  for  rotation  and 

rate-of -strain.  Readers  are  referred  to  the  (1922)  paper  for  the  details  of 

the  ellipsoidal-coordinate  solution.  In  his  equations  (18),  (19)  and  (20)  for 

the  velocity  components,  the  terms  containing  the  constants  A,  B,  C,  F,  F',  G, 

G',  H  and  H'  may  be  rearranged  into  the  stresslet  and  rotlet  distributions  of 

2  2 

equation  (2.7)  plus  octupoles  of  strength  (a^-c  ).  After  some  lengthy 

algebraic  rearrangement0,  these  octupoles  and  the  terms  containing  the 

constants  R,  S,  T,  U,  V  and  W  in  the  (1922)  paper  (these  terms  satisfy 
2 

V  v  =>  0)  reduce  to  the  octupoles  in  equation  (2.7). 


These  steps  are  omitted  here  but  are  available  from  the  author. 


Appendix  2.  COLLECTION  OF  RESULTS  FOR  ELLIPSOIDS  OF  REVOLUTION 


The  scalar  coefficients  which  arise  in  the  solution  of  the  resistance  problems 
for  ellipsoids  of  revolution  are  scattered  throughout  the  literature.  Here, 
the  complete  set  of  resistance  coefficients  for  both  prolate  and  oblate 
spheroids  are  furnished  for  the  convenience  of  the  reader  in  Tables  1-3.  The 
information  is  grouped  as  follows: 

1)  The  expressions  for  the  ellipsoidal  constants,  ou,  3q,  ...  Y£  in  the  limit 
of  oblate  and  prolate  spheroids,  are  given  in  Table  T. 

2)  The  definitions  of  and  expressions  for  the  eight  resistance  functions  which 
relate  the  force,  torque  and  stresslet  on  oblate  and  prolate  spheroids  to 
the  net  translation,  net  rotation  and  ambient  rate-of-strain  are  given  in 
Table  2.  Following  Chwang  &  Wu  (1975),  the  shape-dependence  is  expressed 

in  terms  of  the  eccentricity,  e,  of  the  generating  ellipse. 

3)  In  Table  3,  asymptotic  formulae  are  given  for  all  eight  functions  in  the 
limit  as  e  — >  0  (near-spheres)  and  e  — >  1  (flat  disks  and  thin  needles). 

The  notation  for  the  resistance  functions  follows  that  used  by  Jeffrey  & 

Onishi  (1984).  The  letters  X,  Y  and  Z  are  assigned  according  to  m  -  0,  1  and 

2  respectively,  where  m  is  the  azimuthal  constant  which  appears  in  boundary 

condition.  Superscripts  A,  C,  H  and  M  indicate  the  relation  with  the 

appropriate  resistance  tensor.  The  form  taken  by  tensors  A  and  C  in  Table  2 

is  simply  the  decomposi tion  of  the  translation  and  rotation  problems  into 

motions  parallel  and  perpendicular  to  the  axis  of  symmetry.  The  form  taken  by 

M  is  also  a  consequence  of  the  particle  symmetry.  Finally,  as  a  consequence 

U 

of  the  Lorentz  reciprocal  theorem,  Y  appears  both  as  the  torque  on  a  spheroid 
in  a  rate-of-strain  field  and  also  as  the  stresslet  on  a  rotating  sphere 'd  as 
shown  by  Hinch  (1972). 

The  exact  and  asymptotic  formulae  for  the  resistance  functions  are 
plotted  in  Figures  (4a)  to  (5h),  from  which  it  is  apparent  that  the  asy.qtctic 
forms  are  accurate  over  a  wide  range  of  aspect  ratios.  Rather  curiously,  for 

~  c  u 

a  flat  disk,  all  three  torque  functions,  X",  Y  and  Y“  are  equal  (and  nonzero; 


Oblate  Spheroids 


Ellipsoidal  constants  in  the  limit  of  prolate  and  oblate  spheroidal  cases 
-  y  for  prolate  spheroids  and  B  =  a  for  oblate  spheroids). 
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Table  2  (continued) 

Resistance  Functions  for  Oblate  Spheroids 
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Table  2  (continued) 

Resistance  Functions  for  Prolate  Spheroids 
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Table  2.  The  resistance  functions  for  oblate  and  prolate  spheroids 
(scaled  to  the  spherical  result),  as  a  function  of  e,  the 
eccentricity  of  the  generating  ellipse.  The  constants 

<*l(e),  <>2(e),  Y(e)  and  Y*(e)  are  as  in  Chwang  &  Wu  (1975), 
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Table  3 


Asymptotic  Behavior  of  Oblate  Functions 
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Table  3  (continued) 

Asymptotic  Behavior  of  Prolate  Functions 
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Table  3.  Asymptotic  behavior  of  the  spheroidal  resistance  functions  in 
the  limit  of  spheres,  needles  and  disks. 
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Figure  4.  (Continued) 

e. 

Torque/rate-of-strain  function  Y  . 

f. 

S tress let/rate-of-s train 

(axisyrranetric  straining)  function 

g- 

S tress let/rate-of -strain 

M 

(hyperbolic,  straining)  function  Y 

h. 

Stresslet/ra  te-of-s train 

V* 

(hyperbolic  straining)  function  Z“ 

Aspect  Ratio 


|»A»J 


Aspect  Ratio 


t^ure  5.  The  resistance  functions  for  prolate  spheroids,  scaled  by  the 
results  for  spheres.  The  dashed  curves  are  the  asymptotic 
forms  of  Table  3. 

a.  Force/ translation  (parallel  to  axis)  function  XA. 

b.  Force/ translation  (perpendicular  to  axis)  function  YA. 

p 

c.  Torque/rotation  (parallel  to  axis)  function  X  . 

p 

d.  Torque/rotation  (perpendicular  to  axis)  function  Y  . 


Aspect  Ratio 


Aspect  Ratio 


(Continued) 

e.  Torque/rate-of-strain  function  Y  . 

f.  Stresslet/rate-of-strain  (axisynmetric  straining)  function  X 

M 

g.  Stresslet/rate-of-strain  (hyperbolic  straining)  function  Y  . 

M 

h.  Stresslet/rate-of-strain  (hyperbolic  straining)  function  Z  . 
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that  C  becomes  Isotropic  and 

HmIJ(C"',mnHnkt  reduces  to  **aVd[jkl 
(3.9),  the  Faxen  law  for  the  stresslet  on  a  torque-free  spheroid. 
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